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Flux- averaged and flux-integrated cross sections for quasi-elastic neutrino charged- 
current scattering on nucleus are analyzed. It is shown that the flux-integrated 
differential cross sections are nuclear model-independent. We calculate these cross 
sections using the relativistic distorted-wave impulse approximation and relativistic 
Fermi gas model with the Booster Neutrino Beamline flux and compare results with 
the recent MiniBooNE experiment data. Within these models an axial mass Ma is 
■ extracted from a fit of the measured da/dQ 2 cross section. The extracted value of 

in 

ON , Ma is consistent with the MiniBooNE result. The measured and calculated double 

in 

differential cross sections da / dTd cos generally agree within the error of the exper- 
iment. But the Fermi gas model predictions are completely off of the data in the 



o 



region of low muon energies and scattering angles. 



^ . PACS numbers: 25.30.-c, 25.30.Bf, 25.30.Pt, 13.15.+g 

W. 

I. INTRODUCTION 

The current [1-4] and planed [5] set of accelerator-based neutrino experiments use ex- 
tremely intense neutrino beamlines for precise measurements of the observed neutrino mass 
splitting and mixing angles and detailed experimental study of the neutrino mixing matrix. 
The data of these experiments will greatly extend the statistics due to extremely intense 
neutrino beamline. 

In this situation, the statistical uncertainties should be negligible as compared to sys- 
tematic errors. An important source of systematic uncertainties is related to the neutrino- 
nucleus (vA) cross sections. The neutrino beams of high intensity cover the few-GeV energy 
range. In this energy regime, the dominant contribution to vA cross section comes from 
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charged- current (CC) quasielastic (QE) scattering and resonance production processes. In 
the long-base line neutrino oscillation experiments near and far detectors are used to nor- 
malize the neutrino flux at production and to search for the neutrino oscillation effects. 
While many unknown quantities are eliminated in these experiments by considering ratios 
of far to near events, the cancellation is not complete due to differences in neutrino flux and 
backgrounds in the near and far detectors. Thus, in order to permit precision oscillation 
measurements it is important to have an accurate characterization of the CCQE differential 
cross sections over wide span neutrino energies. 

The current data on CCQE scattering come from a variety of experiment operating at 
differing energies and with different nuclei. The existing data on (anti)neutrino CCQE scat- 
tering come mostly from bubble chamber experiments, which suffer from small statistic. In 
general, the experimental execution and data interpretation are non-trivial for several rea- 
son. Neutrino beams typically span a wide energy range. The neutrino flux itself is often 
poorly known and a background from resonance processes are frequently significant and 
it is difficult for separating from the CCQE signal. Therefore the total QE cross sections 
measured in different experiments with accuracy of 20-40% and even within such large un- 
certainties some results contradict each other. The difference between the total quasielastic 
cross sections, calculated within the framework of various models [8-20] is lower than the 
spread in data. 

More information about neutrino-nuclear CCQE interaction can be obtained from the 
analysis of the charged- current QE event distributions and da/dQ 2 differential cross sections 
as functions of Q 2 (squared four-momentum transfer) [21]. The shape of these distributions 
is sensitive to the Q 2 dependence of two vector, Fi^iQ 2 ), one axial-vector Fa(Q 2 ) form 
factors and nuclear effects. The vector form factors are well-known from electron scattering. 
For the axial-vector form factors the dipole parametrization with one free parameter M A 
(axial mass) is mainly used. This parameter controls the Q 2 dependence of F A (Q 2 ), and 
ultimately, the normalization of the predicted cross sections. The dipole parametrization 
has no strict theoretical basis and the choice of this parametrization is made by the analogy 
with electroproduction. To describe the nuclear effects, neutrino CCQE models typically 
employ a relativistic Fermi gas model (RFGM) [8] in which the nucleons with a flat nucleon 
momentum distribution up to the same Fermi momentum p F and nuclear binding energy e b . 
The experimental values of M A extract from the (anti)neutrino CCQE scattering data, i.e. 
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from the analysis of the shape of the Q 2 - distributions and from the direct measurements of 
the total cross sections. They show very wide spread from roughly 0.7 to 1.2 GeV and the 
resulting world-average M A = 1.03 ± 0.02 GeV [22]. 

Several experiments have recently reported new results on CCQE scattering from high- 
statistics data samples with intense, well-understood neutrino beams. The NOMAD exper- 
iment [23] observe an M A value and cross section (from data taken on carbon) consistent 
with prior world-average. However, data of [1, 24-26] and [27] (preliminary result), collected 
on carbon, oxygen, and iron targets, have indicated a somewhat larger value for M A (by 

10 — 30%). In these experiments the shape of the (^-distribution was analyzed. 

This data show a disagreement with the RFGM predictions. The data samples exhibit 
deficit in the region of low Q 2 < 0.2 (GeV/c) 2 (so-called low- Q 2 problem). As it is known the 
comparison with the low-energy QE electron-nucleus scattering data, the RFGM description 
of this region is not accurate enough [28]. In the region of high-Q 2 the data excess is observed, 
and value of M A) obtained from a fit to the measured data, is higher than the results of 
the previous experiments. The collection of existing results remains puzzling. The next 
experiments MINERvA [6] and MicroBooNE [7] as well as T2K [3] and NOvA [5] near 
detectors will be able to make more precise measurements of the CCQE cross sections in a 
wide range of energies and for various nuclear targets. 

The uncertainties in the theoretical description of the quasielastic neutrino-nucleus scat- 
tering could be considerably reduced if new model-independent absolute differential cross 
section could be provided. The first measurement of the flux-integrated double-differential 
cross section (in muon energy and angle) for CCQE scattering on carbon has been pro- 
duced in MiniBooNE experiment [29]. This cross section contains the most complete and 
model-independent information that is available from MiniBooNE for the CCQE process. 

The aim of this work is to test the RFGM and relativistic distorted-wave impulse ap- 
proximation (RDWIA) predictions against the MiniBooNE data [29]. In the framework 
of these approaches we extract the values of axial mass from the measured flux-integrated 
da/dQ 2 cross section. Then, we calculate with extracted values of Ma the flux- integrated 
differential and flux-unfolded total cross sections and compare the results with data. 

The outline of this paper is the following. In Sec. II we present briefly the RDWIA model 
and discuss the flux-averaged and flux-integrated differential cross sections. The results are 
presented in Sec. III. Our conclusions are summarized in Sec. IV. 
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II. MODEL, FLUX-AVERAGED AND FLUX-INTEGRATED DIFFERENTIAL 

CROSS SECTIONS 

We consider neutrino charged- current QE exclusive 

v{h) + A(p A ) -> nikf) + N{p x ) + B(p B ), (1) 

and inclusive 

uik,) + A(p A ) ^ fi(k f ) + X (2) 

scattering off nuclei in the one-W-boson exchange approximation. Here hi = (£j, fej) and 
fc/ = (ef,kf) are the initial and final lepton momenta, p A = (s A ,p A ), an d Pb = ( £ b,Pb) 
are the initial and final target momenta, p x = (e x ,p x ) is the ejectile nucleon momentum, 
q = (u, q) is the momentum transfer carried by the virtual W-boson, and Q 2 = —q 2 = q 2 —u 2 
is the W-boson virtually. 

A. Model 

The formalism of charged-current QE exclusive and inclusive reactions is described in [18]. 
All the nuclear structure information and final state interaction effects (FSI) are contained 
in the weak CC nuclear tensors W^ v , which are given by bilinear product of the transition 
matrix elements of the nuclear CC operator J M between the initial nucleus state \A) and the 
final state \Bf) as 

WW = T,< B f>pM A )( A \ J t\ B f>P*)> ^ 
f 

where the sum is taken over undetected states. 

We describe CCQE neutrino- nuclear scattering in the impulse approximation (IA), as- 
suming that the incoming neutrino interacts with only one nucleon, which is subsequently 
emitted, while the remaining (A-l) nucleons in the target are spectators. The nuclear cur- 
rent is written as the sum of single-nucleon currents. Then, the nuclear matrix element in 
Eq.(3) takes the form 

(p,S|J"|A) = [ d 3 r exp(it • r)* H (p, r)r"$(r), (4) 
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where T M is the vertex function, t = e B qjW is the recoil-corrected momentum transfer, 
W = a/ {rriA + — Q 2 is the invariant mass, <£> and are relativistic bound-state and 
outgoing wave functions. 

The single- nucleon charged current has V —A structure J M = Jy + J^. For a free-nucleon 
vertex function T M = T v + we use the CC2 vector current vertex function T v . The weak 
vector form factors are related to the corresponding electromagnetic ones for protons and 
neutrons by the hypothesis of the conserved vector current. We use the approximation of 
Ref. [31] on the nucleon form factors. Because the bound nucleons are off shell we employ 
the de Forest prescription [32] and Coulomb gauge for off-shell vector current vertex T v . 
The vector-axial and pseudoscalar form factors are parametrized as a dipole with the axial- 
vector mass, which controls the Q 2 dependence of Fa, and ultimately, the normalization of 
the predicted cross section. 

According to the JLab data [33, 34] the occupancy of the independent particle shell-model 
(IPSM) orbitals of 12 C equals on average 89%. In this work we assume that the missing 
strength (11%) can be attributed to the short-range nucleon- nucleon (NN) correlations in 
the ground state, leading to the appearance of the high-momentum (HM) and high-energy 
component in the nucleon distribution in the target. To estimate this effect in the inclusive 
cross sections we consider a phenomenological model which incorporates both the single- 
particle nature of the nucleon spectrum at low energy (IPSM orbitals) and the high-energy 
and high-momentum components due to NN correlations. 

In the independent particle shall model the relativistic wave functions of the bound nu- 
cleon states $ are obtained as the self-consistent (Hartree- Bogoliubov) solutions of a Dirac 
equation, derived, within a relativistic mean field approach, from Lagrangian containing 
a,u, and p mesons (the a — u model) [35]. We use the nucleon bound-state functions cal- 
culated for carbon by the TIMORA code [36] with the normalization factors S(a) relative 
to the full occupancy of the IPSM orbitals of 12 C: S(lp 3/2 )=84%, S(ls 1/2 )=100%, and an 
average factor of about 89%. These estimations of the depletion of hole states follow from 
the RDWIA analysis of 12 C(e, e'p) for Q 2 < 2 (GeV/c) 2 [34] and are consistent with a di- 
rect measurement of the spectral function using 12 C(e, e'p) in parallel kinematics [39], which 
observed approximately 0.6 protons in a region attributable to a single-nucleon knockout 
fromcorrelated cluster. 

For the outgoing nucleon the simplest choice is to use plane- wave function \1>, i.e., no 
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interactions is considered between the ejected nucleon N and the residual nucleus B (PWIA 
- plane- wave impulse approximation). For a more realistic description, FSI effects should be 
taken into account. In the RDWIA the distorted- wave function \l> are evaluated as solution 
of a Dirac equation containing a phenomenological relativistic optical potential. The channel 
coupling in the FSI [37] of the N + B system is taken into account. The relativistic optical 
potential consists of a real part which describes the rescattering of the ejected nucleon and 
of an imaginary part that accounts for absorption of it into unobserved channels. 

Using the direct Pauli reduction method the system of two coupled first-order radial 
Dirac equations can be reduced to a single second-order Schrodinger-like equation for the 
upper component of Dirac wave function We use the LEA program [38] for the numeri- 
cal calculation of thedistorted wave functions with the EDAD1 parametrization [40] of the 
relativistic optical potential for carbon. This code, initially designed for computing exclu- 
sive proton- nucleus and electron-nucleus scattering, was successfully tested against A(e, e'p) 
data [33, 41] and we adopted this program for neutrino reactions. 

A complex optical potential with a nonzero imaginary part generally produces an ab- 
sorption of the flux. For the exclusive A(/, I'N) channel this reflects the coupling between 
different open reaction channels. However, for the inclusive reaction, the total flux must 
conserve. In Ref. [11, 43] it was shown that the inclusive CCQE neutrino cross section of 
the exclusive channel A(/, I'N), calculated with only the real part of the optical potential is 
almost identical to those of the Green's function approach [11, 42] in which the FSI effects 
on inclusive reaction A(Z, I'X) is treated by means of a complex potential and the total flux 
is conserved. We calculate the inclusive and total cross sections with the EDAD1 relativistic 
optical potential in which only the real part is included. 

The inclusive cross sections with the FSI effects in the presence of the short-range NN 
correlations were calculated using the method proposed in Ref. [18]. In this approach the 
contribution of the NN correlated pairs is evaluated in the PWIA model. We use the 
general expression for the high-momentum and high-energy part of the spectral function 
from Ref. [44] with the parametrization for the nucleon high-momentum distribution from 
Ref. [45], which was renormalized to value of 11%. The FSI effects for the high-momentum 
component is estimated by scaling the PWIA cross section with A(efQf) function determined 
in Ref. [18]. 



7 



B. Flux-averaged and flux-integrated differential cross sections 

In neutrino experiments the differential cross sections of CCQE neutrino-nucleus scatter- 
ing are measured within rather wide ranges of the (anti) neutrino energy spectrum. Therefore 
flux-averaged and flux-integrated differential cross sections can be extracted. 

Because the v^- mode of beams incorporates and spectra the flux-averaged double 
differential cross section a jdTd cos 9 in muon kinetic energy T and muon scattering angle 9 
is the sum of neutrino and antineutrino cross sections 



/ d 2 a \ I d 2 a v \ I d 2 



a 



where 



\dTdcos9 / " \dTdcos9 / \dTdcos9 



\l™ »{T,cos9)) = w v , s (T,cose,ei)—- -(T, cos#, ei)de i: (6) 

\ did cost) I J £l did cos 9 

and W Vi y are weight functions. The normalization of these functions is given by 

f 2 [W U (T, cos 9, £j) + W P (T, cos 9, £j)]d£j = 1. (7) 
The weight functions are defined as follows 

W v , p (T,cos9,e t ) = I uP {ei)/^{T, cos 9), (8) 
where I„ t p is the neutrino (antineutrino) spectrum in u-mode of the flux and 

$(T,cos#) = / 2 [IM) + h(e)]de t (9) 

J ex 

is the neutrino and antineutrino flux which give the contribution to the measured double 
differential cross section at the fixed values of (T, cos#). This flux depends on (T, cos#) 
due to the limits of integration in Eqs. (6), (7) and (9) which are functions of (T, cos#), i.e. 
£i = £mm(T,cos9) and e 2 = e ms , x (T,cos9). In Fig. 1 the double differential cross sections, 
calculates within the RDWIA and RFGM (with the Fermi momentum pp = 221 MeV/c 
and a binding energy = 25 MeV for carbon) are shown as functions of neutrino energy. 
Apparently the ranges [e max (T, cos9) — e min (T,cos9)} where da 2 / dTdcos9 is not equals to 
zero are different in the RDWIA and RFGM. Therefore the value of $(T, cos#) is model- 
dependent and ultimately the weight functions and the cross section (d 2 a /dTd cos 9) depend 
on nuclear models too. Note that the flux $(T, cos 9) should be used to extract the measured 



+ l^ZTn » (5) 
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flux-averaged cross section in the i,j-bms of (T, cos 9) variables (for example see Eq. (3) in 
Ref. [29]). 

Similarly, the flux- averaged da/dQ 2 cross section can be written as sum 

da \ I da v \ I da u \ 

+ ( -17^. ) > ( 10 ) 



dQ 2 / \dQ 2 / \dQ 
where 

/ da u ' p \ f 62 da"' p 

\lQ^ {Q2 ' Tth) ) = J ^(^)^(^ T ^) & - ( n ) 
and T t h is the muon threshold energy after all cuts for CCQE events selection. The weight 
functions in Eq. (11) are defined as follows 

W^,p(QV,) = /„, p(e,)/$(Q 2 ), (12) 

where 

HQ 2 ) = r[Iu(e t ) + U(e)}de t (13) 

is the neutrino and antineutrino flux which gives the contribution to the measured cross 
section at the fixed value of Q 2 . The flux is a function of Q 2 because S\ = e min (Q 2 ) and 
£2 = £max, where e max is the maximal energy in the (anti) neutrino spectrum. The limit 
£min(Q 2 ), and ultimately the flux $(Q 2 ) depend on nuclear model. As a result the extracted 
flux-averaged cross section (da/dQ 2 ) is model-dependent too. 
In Eq.(ll) the cross section da/dQ 2 is defined as 

w (Q2 ' Tth ' £i) ^ FS' w^ (Q2 '^' (14) 

where u cut = min{u ma x{Q 2 ),£i -m^- T th }, m M is the muon mass, u max (Q 2 ) and u min (Q 2 ) 
are the limits of the kinematic allowed w-range at the fixed value of Q 2 . If T th = the upper 
limit cu cut = oo max (Q 2 ). So, the flux-averaged differential (da 2 /dTd cos 9) and (da/dQ 2 ) cross 
sections are model-dependent. 

In Ref. [29] the differential cross sections were extracted using the flux <&bnb that was 
determined by integration the Booster Neutrino Beamline flux [46] over < Si < 3 GeV, i.e. 
$bjvb is a single number (2.90 x 10 11 u^/cm 2 ). Therefore, these flux-integrated differential 
cross sections are not model-dependent and can be written as follows 

/ d 2 a \ mt = t d 2 a» \ int ( d 2 a p V nt 

\dTdcos9j ~ \dTdcos9J + \dTdcos9J ' ^ ' 
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FIG. 1: (Color online) Double differential cross sections vs the neutrino energy calculated in the 
RDWIA (solid line) and RFGM (dashed line) approaches for the four values of (T, cos 9): (0.4 GeV, 
-0.5), (0.4 GeV, 0.7), (0.6 GeV, 0.8), and (1 GeV, 0.8). 

where 

/ d 2 a u ' p \ int T 2 ~ d 2 a u > p 

[iTI^e^V =7 n w^T, a .e,s,) m ^ j (T,cc,e,s t )is l , (ie) 

and 

m) = m) Am • (17) 
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FIG. 2: (Color online) Flux- aver aged (solid line) and flux-integrated (dashed line) double differen- 
tial cross sections versus cos# for T = 0.4 GeV (upper panel) and versus T for cos# = 0.7 (lower 
panel) calculated in the RDWA approach for v- mode of the BNB flux. 

where 

\W {Q2 ' Tth) ) = J w ",'&> e *) s ^® 2 > T *> e *) de '- ( 18 ) 

The weight functions W V} p are defined as 

W v , p (T,COs9,Si) = I v ,u{Si)/^BNB, (19) 
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and 




[W„(T, cos 9, Ei) + W v (T,cos9,e i )]de i < 1. 



(20) 



because of $bnb > $(T, cos#) and §bnb > &(Q 2 )- These functions depend only on 
(anti) neutrino energy and are model-independent. As an example, on Fig. 2 the flux- 
averaged and flux-integrated double differential cross sections calculated within the RDWIA 
model for v- mode of the BNB flux are compared. Apparently the flux-averaged cross sec- 
tions are higher than the flux-integrated ones. This is because the normalization of W Uj p 
(Eq.(20)) is less than unit. From the practical point of view, the flux- integrated differential 
cross sections are more useful than flux-averaged ones because they are not model-dependent 
and can be used for comparison to models of CCQE interaction on nuclear targets. 



New data for CCQE events (^-distribution measured in the MiniBooNE experiment 
were presented in Refs. [29, 30]. The CC one pion production (CCl7r + ) background was 
measured and subtracted instead of calculated one [1]. With measured CCl7r + background 
incorporated, a "shape-only" fit to the CCQE events sample was performed to extract values 
for adjusted CCQE model parameters, M A and k within the Fermi gas model. To tune this 
model to the low Q 2 , the parameter k was introduced [1] which reduced the phase space 
volume at low-momentum transfer. Note that at k — 1 the phase space volume is the same 
as well as in the "standard" RFGM. This parameter controls the Q 2 - distribution only in the 
low-Q 2 region. The shape-only fit yields the model parameters, M A = 1.35 ± 0.17GeV/c 2 
and k = 1.007 ± 0.012. The extracted value for M A is approximately 30% higher than the 
world averaged one. 

The MiniBooNE CC flux- integrated single differential cross section da/dQ 2 per neu- 
tron was extracted as a function of Q 2 in the range < Q 2 < 2 (GeV/c) 2 . To extract value 
for the parameter M A we calculated this cross section with the BNB flux in the RDWIA 
and RFGM models using the Q 2 -bins AQ 2 = Q 2 i+1 - Q 2 similar to Ref. [29] 



III. RESULTS AND ANALYSIS 



A. CCQE flux-integrated da/dQ 



2 differential cross section 




(21) 
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FIG. 3: (Color online) Flux-integrated da/dQ 2 cross section per neutron target for the CCQE 
scattering. Calculations from the RDWIA with M A = 1.37 GeV/c 2 and RFGM with M A = 1.36 
GeV/c 2 . The MiniBooNE data are shown as points with the shape error only. 

Because the data include events with T M < 200 MeV [29], we calculated da/dQ 2 with 
7) t = in Eq. (14). Within the RDWIA (RFGM) model the fit to the extracted flux- 
integrated da/dQ 2 cross section with only-shape error yields the parameter M a = 1.37±0.05 
GeV/c 2 (M A = 1.36±0.05 GeV/c 2 ). These values are consistent with the MiniBooNE result 
M A = 1.37 ±0.17 GeV/c 2 . 
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Figure 3 shows measured flux- integrated da/dQ 2 differential cross section as a function 
of Q 2 compared with the RDWIA (M A = 1.37 GeV/c 2 ) and RFGM (M A = 1.36 GeV/c 2 ) 
calculations. There is an overall agreement between the RDWIA result and the data across 
of the full range Q 2 = ± 1 (GeV/c) 2 , whereas the RFGM overestimates the measured 
differential cross section at Q 2 < 0.2 (GeV/c) 2 . At higher Q 2 a good match between the 
RFGM calculated and measured cross sections is observed. Thus, so-called low-Q 2 problem 
is successfully solved in the distorted-wave approach. 



B. CCQE flux-integrated double differential cross section 

The flux-integrated double differential cross section per neutron d 2 a/dTdcos9, for the 
CCQE process was extracted in Ref. [29] for the kinematic range, —1 < cos9 < 1, 
0.2 < T < 2 GeV. The flux-integrated CCQE total cross section, obtained by integrating 
the double differential one over this range was measured to be 18.447 x 10~ 39 cm 2 and 
9.429 x 10~ 39 cm 2 for range -1 < cos 9 < 1, < T < 2 GeV. The total normalization 
error on this measurement is 10.7%. These results contain the most complete and model- 
independent information that is available from experiment on the CCQE process. 

We calculated the flux-integrated double differential cross section (d 2 a /dTd cos 6) int ) for 
the BNB v mu flux within the RDWIA and RFGM models with the extracted values of Ma 
using the T and Q 2 -bins similar to Ref. [29] 



(fa \ int l fTi+i /-(cos 9)^+1 



dTd cos 9 J ATA cos 9 JTi , (cos e) . 



d 2 a 



dTd cos 9 



(T, cosfl) 



int 



dTd cos 9, (22) 



where AT = T i+1 -T t = 0.1 GeV and A cos 9 = (cos9) j+1 - ( y cos9) j = 0.1. 

Figures 4 and 5 show measured flux-integrated d 2 a/dTdcos9 cross sections as functions 
of cos6> for several bins of muon kinetic energy in the range 0.2 < T < 2 GeV as compared 
with the RDWIA and RFGM calculations. There is good agreement between the RDWIA 
calculations and data within the error of the experiment. But in the regions 0.2 < T < 0.3 
GeV, -1 < cos# < -0.3 and 0.2 < T < 0.5 GeV, 0.9 < cos# < 1 the RDWIA result are 
slightly lower then measured cross section and the difference decreases with muon energy. 

The RFGM prediction also agree well with data within the errors, except the region 
0.7 < cos# < 1 and 0.2 < T < 0.5 GeV, where the calculated cross sections fall down 
rapidly with cos#. In this kinematic region the Fermi gas model underestimates the double 
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FIG. 4: (Color online) Flux-integrated d 2 a /dTdcos 6 cross section per neutron target for the za, 
CCQE process as a function of cos# for the four muon kinetic energy bins: T(GeV)=(0.2 - 0.3), 
(0.3 - 0.4), (0.4 - 0.5), and (0.5 - 0.6). As shown in the key, cross sections were calculated within 
the RDWIA (M A = 1.37 GeV/c 2 ) and RFGM (M A = 1.36 GeV/c 2 ). The MiniBooNE data are 
shown as points with the shape error only. 
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FIG. 5: (Color online) Same as Fig. 4 but for muon kinetic energy bins: T(GeV)=(0.6 - 0.7), (0.7 
- 0.8), (0.8 - 0.9), and (0.9 - 1). 

differential cross section significantly. This trend is characteristic of nucleon momentum 
distribution and Pauli bloking effect as calculated in the Fermi gas model [47]. 

Fifure 6 shows measured flux- integrated d 2 a jdTd cos 9 cross sections as functions of muon 
energy for four bins of muon scattering angle as compared with the RDWIA and RFGM 



16 




FIG. 6: (Color online) Flux-integrated d 2 a /dTdcos 6 cross section per neutron target for the za, 
CCQE process as a function of muon energy for the four muon scattering angle bins: cos #=(0.6 
- 0.7), (0.7 - 0.8), (0.8 - 0.9), and (0.9 - 1). As shown in the key, cross sections were calculated 
within the RDWIA (M A = 1.37 GeV/c 2 ) and RFGM (M A = 1.36 GeV/c 2 ). The MiniBooNE data 
are shown as points with the shape error only. 
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calculations. Apparently that the RDWIA cross sections are lower than the measured ones 
in the kinematic region 0.9 < cos6> < 1, 0.2 < T < 0.5 (GeV) and the RFGM calculation 
underestimates the measured double differential cross section significantly in the range 0.7 < 
cos 9 < 1, 0.2 < T < 0.5 GeV. 

So, the comparison measured and calculated flux-integrated d 2 a /dTd cos 9 cross sections 
shows that the Fermi gas model prediction are completely off the data in the range 0.7 < 
cos 6* < 1, 0.2 < T < 0.5 GeV. The RDWIA cross sections underestimate the measured ones 
for muon production with energies T < 0.3 GeV and scattering angles cos^ > 0.9. 

C. CCQE flux-integrated da/dT and da/dcosd cross section 

The flux-integrated single differential cross sections da/dT and da/dcos9 (for T > 0.2 
GeV) are presented in Fig. 7, which shows da/dT as a function of kinetic muon energy and 
da/d cos 9 as a function of muon scattering angle. Here the results obtained in the RDWIA 
and Fermi gas models compared with the MiniBooNE data. The measured flux-integrated 
da/dT (da/d cos 9) cross section with the shape error has been obtained by summing the 
double differential one over cos#-bins (T-bins) presented in Tables VI and VII in Ref. [29]. 
There is a good agreement between the calculated and measured cross sections, with the 
excepton of the bin 0.2 < T < 0.3 GeV. The flux integrated total cross sections obtained 
in the RDWIA and RFGM approaches by integrating the double differential cross sections 
(over -1 < cos0 < 1, 0.2 < T < 2 GeV), are equal of 8.208 x 10~ 39 cm 2 and 8.310 x 10~ 39 
cm 2 , correspondingly, and agree with measured one of 8.447 x 10~ 39 cm 2 . 

D. CCQE total cross section 

The MiniBooNE flux-unfolded CCQE cross section per neutron as a function of neutrino 
energy is shown in Fig. 8 together with data of Refs. [48-51]. Also shown for comparison are 
the results obtained in the RDWIA, PWIA, and RFGM approaches. The calculated cross 
sections with the values of M A) extracted from the shape-only fit to the flux-integrated da/Q 2 
data reproduce the MiniBooNE total cross section within the errors of experiment over the 
entire measured energy range. At the average energy of the MiniBooNE flux 800 MeV), 
the extracted cross section is ~ 30% higher than what is commonly assumed for this process 
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FIG. 7: (Color online) Flux-integrated da/dT cross section as a function of muon energy (upper 
panel) and da jd cos 6 cross section for T > 0.2 GeV as a function of muon scattering angle (lower 
panel) for the i/„ CCQE process. As shown in the key, cross sections were calculated within the 
RDWIA and RFGM. The MiniBooNE data are shown as points with the shape error only. 
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FIG. 8: (Color online) Total v„ CCQE cross section per neutron as a function of neutrino energy. 
Data points for different targets are from [26, 29, 48-51]. Also shown are predictions of the RDWIA 
(M A = 1.37 GeV/c 2 ), PWIA (M A = 1.37 GeV/c 2 ), and RFGM (M A = 1.36 GeV/c 2 ). 

assuming the RFGM and world-average value of the axial mass, M A = 1-03 GeV/c 2 . Note, 
that the spread in the data is much higher than a difference in predictions of the RDWIA, 
PWIA, and RFGM approaches. So, the comparison of the predicted and measured model- 



20 



independent flux-integrated double differential cross sections is more sensitive test of the 
employed models of the CCQE process than the comparison of the total cross sections. 

IV. CONCLUSIONS 

In this paper, we analyze the flux-averaged and flux-integrated differential and total 
CCQE cross sections placing particular emphasis on their nuclear-model dependence. We 
found that the flux-integrated cross sections are model-independent and can be used to test 
of employed models of the CCQE interaction on nuclear targets. The flux-integrated double 
differential d 2 a / dT d cos 0, single differential da/dQ 2 , da/dT, da/dcos9, and flux-unfolded 
a{ei) CCQE cross sections were measured in the MiniBooNE experiment [29]. 

Using the RDWIA and RFGM approaches with the BNB flux we extracted an axial 
mass from a "shape-only" fit of the measured flux-integrated da/dQ 2 differential cross sec- 
tion. The extracted value of M A = 1.37 ± 0.05 GeV/c 2 (RDWIA) and M A = 1.36 ± 0.05 
GeV/c 2 (RFGM) that is consistent with the MiniBooNE result of Ma = 1.35±0.17 GeV/c 2 . 
The flux-integrated double differential cross sections were calculated in these models with 
extracted values of Ma- There is an overall agreement between the RDWIA result and 
data, whereas the RFGM calculation overestimates the measured cross section at Q 2 < 0.2 
(GeV/c) 2 . Thus, so-called low-Q 2 problem is successfully solved in the framework of RD- 
WIA. 

We also calculated in the RDWIA and RFGM approaches the flux-integrated 
d 2 a/dTdcos9, da/dQ 2 , da/dT (for muons with kinetic energy T > 0.2 GeV), and total 
cross sections and compared them with the MiniBooNE data. The comparison of the RD- 
WIA double differential cross section shows good agreement with data within the error of 
the experiment, except in the region 0.2 < T < 0.3 GeV, 0.9 < cos 9 < 1 where the calcu- 
lated cross sections are lower then measured ones. A good agreement between the RFGM 
calculation and data is observed exclusive of the range 0.7 < cos 01, 0.2 < T < 0.5 GeV 
where the Fermi gas model predictions are completely off of the data. The calculated da/dT 
and da /d cos 9 also describe well the measured cross sections except the muon energy bin 
0.2 < T < 0.3 GeV where the calculations are lower then the data. 

The calculated and measured flux-integrated total cross sections are match well. The RD- 
WIA, PWIA and RFGM calculations with extracted values of Ma reproduce the MiniBooNE 
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flux-unfolded CCQE cross section within the experimental error over the entire measured 
energy range. 

We conclude that the flux-integrated double differential cross section which is model- 
independent should be used as the preferred choice for comparison to employed model of 
the CCQE interaction on nuclear targets. 
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